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Abstract 

In this talk we present the results obtained from a study of 0(2000) (quenched) 
lattice configurations from the APE collaboration, at 6.0 < P < 6.4, using both the 
Wilson and the SW-Clover fermion action. We determine the light hadronic spectrum 
and the meson decay constants. For the light-light systems we find an agreement with 
', the experimental data of ~ 5% for mesonic masses and ~ 10% — 15% for baryonic 

masses and pseudoscalar decay constants; a larger deviation is present for the vector 
decay constants. For the heavy-light decay constants we find fo^ = 237 it 16 MeV, 
Jd = 221 ± 17 MeV {fos/fD = 1-07(4)), fs, = 205 ± 35 MeV, Jb = 180 ± 32 MeV 

^ ' Introduction 

I The aim of this exposition is to describe the recent results obtained from the APE col- 

' laboration on the hadronic spectrum and the meson decay constants. A full description 

■ of the method used and the results obtained are reported in [||, ^. 

I The lattice technique has proved a very effective theoretical tool to determine phe- 

f-*) I nomenological quantities such as the mass spectrum and weak decay matrix elements. 

' Unlike other approaches, it does not (in principle) suffer from uncontrolled approxima- 

, tions. However, in practice, one is forced to work on a lattice of (i) finite size, with (ii) 

"j^ I a finite lattice spacing, and (iii) unphysically large masses for the light quarks. Also 

■ the quenched approximation is often used in lattice studies. 

Oh. The aim of the high statistics simulations [|l|, ^ is to study the main systematic errors 

(— I I present in the extraction of the light hadronic spectrum and meson decay constants, 

sjjl ■ It is important to study the systematics due to (i), (ii) and (iii) before the effects of 

•rH , the quenched approximation can be correctly understood. 

r> i The results reported in this talk are obtained from 0(2000) (quenched) lattice config- 

' urations from the APE collaboration, for different lattice volumes at 6.0 < (3 < 6.4 

using both the Wilson action and the SW-Clover fermion action. The main parameters 
used in each simulation are listed in table 1. The values of beta have been chosen: 

a) small enough to obtain accurate results on reasonably large physical volumes; 

b) large enough to be in the scaling region. 

In this range of f3 we cannot draw any conclusion about a dependence of hadron masses 
and pseudoscalar decay constants. On the other hand, we find that the quenched ap- 
proximation gives a resonable agreement with the experimental values, when the com- 
parison is possible. 

The main physical results of our study have been given in the abstract. 
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Table 1: Summary of the parameters of the runs analyzed. 





C60b 


C60a 


W60 


C62a 


W62a C62b 


W62b 


W64 


C64 




6.0 


6.0 


6.0 


6.2 


6.2 6.2 


6.2 


6.4 


6.4 


Action 


sw 


SW 


Wil 


SW 


Wil SW 


Wil 


Wil 


SW 


# Confs 


170 


200 


320 


250 


250 200 


110 


400 


400 


Volume 


18^ X 64 18^ X 64 18^ x 64 24^ x 64 24^ x 64 18^ x 64 24^ x 64 24^ x 64 24^ x 64 


Ki 








0.14144 


0.1510 




D 1488 


D i4nn 




0.1425 


0.1425 


0.1530 


0.14184 


0.1515 0.14144 


0.1510 


0.1492 


0.1403 




0.1432 


0.1432 


0.1540 


0.14224 


0.1520 0.14190 


0.1520 


0.1496 


0.1406 




0.1440 


0.1440 


0.1550 


0.14264 


0.1526 0.14244 


0.1526 


0.1500 


0.1409 




0.1150 




0.1255 


0.1210 


0.1300 


0.1300 








0.1200 




0.1320 


0.1250 


0.1350 


0.1350 








0.1250 




0.1385 


0.1290 


0.1400 


0.1400 








0.1330 




0.1420 


0.1330 


0.1450 


0.1450 












0.1455 






0.1500 










Li| 


jht-light 


mesons with zero momentum 








15-28 


15-28 


15-28 


18-28 


18-28 18-28 


18-28 


24-30 


24-30 






Heavy-light 


mesons 


with zero momentum 






tl-t2 


15-28 


15-28 


15-28 


20-28 


20-28 


20-28 












Baryons with zero momentum 








tl-t2 




12-21 


12-21 


18-28 


18-28 18-28 


18-28 


22-28 


22-28 


"1^ 


2.00(10) 


1.98(8) 


2.26(5) 


2.7(1) 


3.00(9) 3.0(3) 


3.0(1) 


-1.1(2) 


-1.0(2) 



Lattice Details 

Hadron masses and decay constants have been extracted from two-point correlation 
functions in the standard way. For the meson masses and decay constants we have 
computed the following propagators: 

GssW = $^(P5(X,0^'5^(0,0)) , 

X 

Go5(t) = ^(Ao(x,i)P5^(0,0)) , (1) 

X 

where 

P5{x,t) = iq{x,t)-f5q{x,t) , 
A^{x,t) = q{x,t)-ff,-f5q(x,t) . 

and the following propagators of the vector states: 

Gu{t)= E T.{V^{x,t)V^{0,0)), (2) 

1=1,3 X 

where 

Viix) = q{x,t)-/iq{x,t) . 
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In order to determine the baryon masses we have evaluated the following propagators: 



Gn{t) = 5](iV(x,t)iVt(o,o)) , 

X 

Gsit) = ^(A^(x,t)Aj,(0,0)), (3) 



where 



N = eaUy"G^^S')u^ 

We have fitted the zero-momentum correlation functions in eqs. |l|, ^ and |^ to a single 
particle propagator with cosh or sink in the case of mesonic and axial-pseudoscalar 
correlation functions and with an exp function in the case of the baryonic correlation 
functions 

^55 -j^ 'J-' 

GssW = ^exp(--Mp5T)cosh(Mp5(--t)) , 



Gu{t) = — exp(--MyT)cosh(My(--t)) 



^05 

Go5{t) = j^eM-^MpsT)smh{Mps{--t)) 



(4) 



GnAt) = C"'^exp(-M„,5t) , 

in the time intervals reported in table 1. In eqs. ^, T represents the lattice time 
extension, the subscripts PS and V stand for pseudoscalar and vector meson, n and 
6 stand for nucleon- and delta-like baryons^. To improve stability, the meson (axial- 
pseudoscalar) correlation functions have been symmetrized (anti-symmetrized) around 
t = T/2. The time fit intervals have been chosen with the following criteria: we fix the 
lower limit of the interval as the one at which there is a stabilization of the effective 
mass, and, as the upper limit the furthest possible point before the error overwhelms 
the signal. 

The pseudoscalar and vector decay constants fps and 1/fv are defined through the 
equations 

{0\Ao\PS) = i^Mps, (5) 

mi\y.r) = e\-^, (6) 
JV^V 

where e[ is the vector-meson polarization, Mpg and My are the pseudoscalar and 
vector masses and Zy^A are the renormalization constants. We have extracted fps 



^i.e. n stands for the nucleon TV, AS or ^ baryon, and 5 is either a or a VL. 
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Table 2: Predicted hadron masses in GeV for all lattices, using the scale from Mk*- 





Mp 


Ms 


Mn 


Mas 


Mh 


Ma 


Mn 


Exper. 


0.770 


1.019 


0.9389 


1.135 


1.3181 


1.232 


1.6724 


C60a 


0.809(7) 


0.977(7) 


1.09(5) 


1.21(4) 


1.32(4) 


1.3(1) 


1.60(9) 


W60 


0.808(3) 


0.978(3) 


1.19(5) 


1.29(4) 


1.40(4) 


1.46(7) 


1.71(4) 


C62a 


0.81(1) 


0.98(1) 


1.1(1) 


1.22(8) 


1.34(7) 






W62a 


0.803(6) 


0.984(6) 


1.17(7) 


1.28(6) 


1.39(5) 






C62b 


0.79(1) 


1.00(1) 


1.1(2) 


1.2(2) 


1.4(1) 


1.6(3) 


1.9(2) 


W62b 


0.797(7) 


0.989(7) 


1.2(1) 


1.3(1) 


1.40(9) 


1.50(9) 


1.72(5) 


W64 


0.796(4) 


0.990(4) 


1.21(9) 


1.32(8) 


1.43(6) 


1.4(2) 


1.72(9) 


C64 


0.792(4) 


0.994(4) 


1.2(1) 


1.29(8) 


1.41(7) 


1.3(2) 


1.7(1) 



from the ratio 



^05 2^ 

ZA^tanh{Mps{--t)) 

^^^^^^^tanh(Mp5(?-t)) , (7) 



VZ^ ' '2 

and the vector-meson decay constant has been obtained directly from the parameters 
of the fit to Gii{t), eqs. 13: 



1 Vz 



^- (8) 



Zvfv My 

Light-Light hadron masses and decay constants 

In order to reduce the error coming from the chiral extrapolation, we have extracted 
as much physics as possible from the "strange" region. The method we have used to 
extract the hadron masses and the light-light mesons decay constants from the lattice 
data is fully explained in [|l[] . In tables § - we give a list of the lattice predictions 
for the light mass spectrum and light-light decay constants from C60a - C64. The scale 
is set from the K* mass. Note that all errors quoted in tables |2| - ^ are statistical only. 
As can be seen from table for the hadron masses there is good consistency of the 
physical predictions between the different simulations within the errors. 
To compare the lattice decay constants with the experimental ones, we have used a 
'boosted' one-loop form of the renormalization constants @ - 0|: 



Wilson action Za = 1 - O.lSAgj 



Zv = 1 - 0.174C/2 



MS 
2 

MS 
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Table 3: Predicted meson decay constants for all lattices, using the scale from Mk* 



Exper. 



U (GeV) 



0.1307 



0.28 



fK {GeV) 
0.1598 



fx* 



0.23 



C60a 
W60 
C62a 
W62a 
C62b 
W62b 
W64 
C64 



0.134(9) 
0.155(7) 
0.124(9) 
0.135(6) 
0.14(2) 
0.135(8) 
0.147(9) 
0.144(9) 



0.35(3 
0.37(2 
0.33(3 
0.35(2 
0.26(4 
0.36(2 
0.29(1 
0.25(1 



0.149(8) 
0.167(6) 
0.143(8) 
0.153(5) 
0.16(2) 
0.152(7) 
0.161(8) 
0.158(8) 



0.33(2) 
0.35(1) 
0.30(2) 
0.33(1) 
0.25(3) 
0.34(2) 
0.283(9) 
0.25(1) 



0.30(1) 
0.324(8) 
0.281(7) 
0.307(6) 
0.25(2) 
0.315(9) 
0.272(6) 
0.243(8) 



Clover action Za = 1 — 0.0177(7|j^ 
Zy =l-Q.mg^ 

where = 6//%^, with 

/?M5=<^,/ag>/3 + 0.15. (9) 

The results are shown in table ^. For the ratio fps/J^v we see again a global consis- 
tency of our data both for Wilson and the SW-Clover action. We do not see a strong a 
dependence for either actions comparing lattices at different (5 (C60a, C62a and W60, 
W62a). The experimental points lie quite well in the extrapolated/interpolated lattice 
data. 

The values of the pseudoscalar decay constants in table ^ show larger deviations from 
the experimental data than the ratio fps/My- It is clear that a cancellation of system- 
atic error occurs in this ratio. Still there is an agreement within 1.5 standard deviations 
with the experimental data apart from the W60 value which is quite high. 
The situation is more delicate for the vector decay constant for which a dependence 
on the volume and may be present but with our data it would be difficult to dis- 
entangle the two effects [|l|. The 'strange' vector decay constants seem to be slightly 
more stable also because no extrapolation is needed. 



Strange Quark mass 

Lattice QCD is in principle able to predict the mass of a quark from the experimental 
value of the mass of a hadron containing that quark. The 'bare' lattice quark mass m(a) 
can be extr acted directly from lattice simulations and can be related to the continuum 
mass m^^'^(/x) renormalized in the minimal-subtraction dimensional scheme through a 
well-defined perturbative procedure (5|. Following ref. 

m^(^) = z'S^{f,a)m{a) 
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Table 4: Values for the lattice strange lattice quark masses for all lattices and the corre- 
sponding MS values at NLO, both in MeV . 





' ' \ J 


mf^{pL = 2 GeV) 


C60a 


89(3) 


120(10) 


W60 


98(2) 


130(20) 


C62a 


83(4) 


120(10) 


W62a 


93(3) 


130(10) 


C62b 


75(6) 


110(10) 


W62b 


92(4) 


130(20) 


W64 


82(3) 


120(10) 


C64 


69(3) 


100(10) 



where m(a) is the bare lattice quark mass and (//o) is the mass renormalization 
constant at scale ^ which we choose to be 2 GeV . The bar e an d renormalized strange 
quark masses are reported in table ^. The error on m^'^^{pL) has been estimated 
as in ref. ||5| taking into account the spread due to different definitions of the strong 
coupling constant. There is a good consistency among the values coming from the 
different lattices and we do not see any dependence on the lattice spacing a within 
the errors in the Wilson data and a mild tendency in the clover data to decrease with 
increasing fi. We therefore conclude that any 0{a) effects present are beneath the level 
of statistics, and/or hidden among finite volume effects. We then extract the average 
value of the strange quark mass without any extrapolation and get: 

m^i^i = 2GeV) = 122 ± 20MeF 

which is in agreement with the result of ref. ^j. It is also compatible with the value 
mf^(// = 2GeV) = 100 ± 21 ± 10 MeV of ref. [§, but one should take into account 
that this value comes from an analysis on various Wilson and Staggered lattices at 
different values of (3 and an extrapolation in a. 



Heavy-Light decay constants 

The major sources of uncertainty in the determination of the heavy-light pseudoscalar 
decay constants, besides the effects due to the use of the quenched approximation, 
come from the calculation of the constant Za in eq. (^) and from discretization errors 
of 0{a) present in the operator matrix elements. A method to get rid of Za consists of 
extracting the decay constants of heavier pseudoscalar mesons by computing the ratio 
Rp = fpjfT^ and multiplying Rp by the experimental value of the pion decay constant. 
Hopefully, by taking the ratio, some of the 0(a) effects are eliminated. These effects 
are expected to be more important for jp, than for f^^^K since, at current values of a, 
the relevant parameter mcharmO^ is not very small. A full analysis of the discretization 



errors is performed in the section 2 of [M. The results are reported in table 5 . We will 



^For the run W60, the heavy-light correlators are computed on a subset of 120 configurations. 
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Table 5: Predicted heavy-light meson decay constants for all lattices. 



Run 


fos 1 Ik 


Id/ f-K 


fos 1 fn 


Ibs I Ik 




fSs 1 Ib 


C60b 


1.56(3) 


1.63(4) 


1.08(1) 


1.48(7) 


1.53(9) 


1.10(3) 


C62a 


1.48(6) 


1.58(8) 


1.07(2) 


1.28(9) 


1.32(13) 


1.14(6) 


W603 


1.11(3) 


1.14(4) 


1.06(1) 


0.79(4) 


0.81(5) 


1.05(2) 


W62a 


1.23(4) 


1.31(6) 


1.07(1) 


0.83(4) 


0.86(6) 


1.10(3) 


W62b 


1.19(5) 


1.25(7) 


1.09(2) 


0.84(5) 


0.86(6) 


1.12(3) 



first discuss the results for charmed mesons, for which the extrapolation in the heavy 
quark mass is not a relevant source of systematic uncertainties, and then discuss the B- 
meson case. We find that higher statistics and larger intervals of a values are required 
to satisfactorily uncover the o dependence of the decay constants. An extrapolation to 
a = it is not possible at this stage. Thus we believe that the best estimate of fo^ is 
obtained from the Clover data at /? = 6.2, by using Rds = fosl Ik (from a linear fit in 
the light quark masses, a quadratic fit in and without any KLM factor). As for 

the error, we take as a conservative estimate of the discretization error the difference 
between the results obtained at /? = 6.0 and 6.2, and combine it in quadrature with 
the statistical one. This gives Rd^ = 1.48(10) from which, by using f^"^ = 159.8 MeV, 
we obtain fu^ = 237 it 16 MeV. This value is in very good agreement with the ex- 
perimental value 00 of branching fraction B{Df fJ-'^i^fj.) = (4.6 ± 0.8 ± 1.2) x 10~^ 
which corresponds to a decay constant value of fo^ = (241 it 21 it 30) MeV (note 
that was predicted by lattice calculations long before its measurement). By using 
Rd, = 1.48(10) combined with /dJId = 1.07(4) we obtain fo = 221 ± 17 MeV. We 
believe that these are our "best" results. 

In order to obtain fs^ and Jb, an extrapolation in the heavy quark mass well outside 
the range available in our simulations is necessary. Discretization errors can affect the 
final results in two ways. Not only they can change the actual values of the decay 
constants, but also deform the dependence of fp on uih- Moreover, points obtained 
at the largest values of mna become the most important, since we extrapolate in the 
direction of larger values of uih ■ A full analysis of these effects is done in . 
In order to extract our "best" values for fp^ and fp we have proceed exactly as for the 
D-meson case. We have obtained = 205 it 15 ± 31 MeV= 205 it 35 MeV, where the 
second error (31 MeV) is the discretization error, estimated by comparing the results 
from C60b and C62a as done for fo,. We also have obtained /b = 180 ± 32 MeV. The 
decay constants of the B-mesons are not yet measured and the numbers given above 
are predictions of the lattice. The results we have obtained for the heavy-light decay 
constants are in very good agreement with the compilations of lattice calculations pre- 
sented in refs. pi p^. 



"^There is a new measurement of B{D^ r Vr) from L3 They have obtained B{D^ t Vr) = 
(0.074 ± 0.028(stat) ± 0.016(syst) ± 0.018(norm)) which corresponds to fo, = (309 ± 58(stat) ± 33(syst) ± 
38 (norm)) MeV. 
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Conclusions 



In this talk we have reported results from a large set of data on lattices of different 
lattice spacing and lattice sizes obtained with both the Wilson and SW-Clover actions. 
The results obtained and reported in the abstract give us confidence in the ability of 
the lattice to predict other non-perturbative quantities that are of phenomenological 
interest. Further studies, with comparable (or smaller) statistical errors and physical 
volume, at smaller values of the lattice spacing, are required to reduce the systematic 
error due to 0(a) effects. The systematic errors will be completely understood only 
when the quenched approximation will be removed. 
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